Coping with salinity extremes: Gill transcriptome profiling in the black-chinned tilapia (Sarotherodon melanotheron).

Steeper and sometimes extreme salinity gradients increasingly affect aquatic organisms because of climate change. Hypersalinity habitats demand powerful physiological adaptive strategies. Few teleost species have the capacity to spend their whole life cycle in salinities way over seawater levels. Focusing on the multifunctional gill, we unraveled the tilapia S. melanotheron key strategies to cope with different environmental conditions, ranging from freshwater up to hypersaline habitats. De novo transcriptome assembly based on RNAseq allowed for the analysis of 40,967 annotated transcripts among samples collected in three wild populations at 0, 40 and 80 ‰. A trend analysis of the expression patterns revealed responses across the salinity gradient with different gene pathways involved. Genes linked to ion transport, pH regulation and cell surface receptor signaling were mainly upregulated in the high salinity habitat. We identified tight junction proteins that were critical in high salinity habitats and that were different from the well-known tightening junctional proteins identified and expressed in fresh water. Expression profiles also suggest a change in the vascular tone that could be linked to an osmorespiratory compromise not only in fresh water, but also in high salinity environments. A striking downregulation of genes linked to the immune system and to the heat shock response was observed suggesting an energetic trade-off between immunity and acclimation/adaptation in the hypersaline habitat. The high expression of transcripts coding for immune and heat shock response in the freshwater habitat suggests the establishment of powerful mechanisms to protect gills from environmental threats and to maintain protein integrity. Non-directional expression trends were also detected with an upregulation of genes only in the hypersaline habitat (80 ‰) or only in the marine habitat (40 ‰). Unravel physiological strategies in S. melanotheron populations will help to better understand the molecular basis of fish euryhalinity in salinity-contrasted environments.

Genome-wide analysis of European sea bass provides insights into the evolution and functions of single-exon genes.

Several studies have attempted to understand the origin and evolution of single-exon genes (SEGs) in eukaryotic organisms, including fishes, but few have examined the functional and evolutionary relationships between SEGs and multiple-exon gene (MEG) paralogs, in particular the conservation of promoter regions. Given that SEGs originate via the reverse transcription of mRNA from a "parental" MEGs, such comparisons may enable identifying evolutionarily-related SEG/MEG paralogs, which might fulfill equivalent physiological functions. Here, the relationship of SEG proportion with MEG count, gene density, intron count, and chromosome size was assessed for the genome of the European sea bass, Dicentrarchus labrax. Then, SEGs with an MEG parent were identified, and promoter sequences of SEG/MEG paralogs were compared, to identify highly conserved functional motifs. The results revealed a total count of 1,585 (8.3% of total genes) SEGs in the European sea bass genome, which was correlated with MEG count but not with gene density. The significant correlation of SEG content with the number of MEGs suggests that SEGs were continuously and independently generated over evolutionary time following species divergence through retrotranscription events, followed by tandem duplications. Functional annotation showed that the majority of SEGs are functional, as is evident from their expression in RNA-seq data used to support homology-based genome annotation. Differences in 5'UTR and 3'UTR lengths between SEG/MEG paralogs observed in this study may contribute to gene expression divergence between them and therefore lead to the emergence of new SEG functions. The comparison of nonsynonymous to synonymous changes (Ka/Ks) between SEG/MEG parents showed that 74 of them are under positive selection (Ka/Ks > 1; p = .0447). An additional fifteen SEGs with an MEG parent have a common promoter, which implies that they are under the influence of common regulatory networks.

Thermal selection as a driver of marine ecological speciation.

Intraspecific genetic structure in widely distributed marine species often mirrors the boundaries between temperature-defined bioregions. This suggests that the same thermal gradients that maintain distinct species assemblages also drive the evolution of new biodiversity. Ecological speciation scenarios are often invoked to explain such patterns, but the fact that adaptation is usually only identified when phylogenetic splits are already evident makes it impossible to rule out the alternative scenario of allopatric speciation with subsequent adaptation. We integrated large-scale genomic and environmental datasets along one of the world's best-defined marine thermal gradients (the South African coastline) to test the hypothesis that incipient ecological speciation is a result of divergence linked to the thermal environment. We identified temperature-associated gene regions in a coastal fish species that is spatially homogeneous throughout several temperature-defined biogeographic regions based on selectively neutral markers. Based on these gene regions, the species is divided into geographically distinct regional populations. Importantly, the ranges of these populations are delimited by the same ecological boundaries that define distinct infraspecific genetic lineages in co-distributed marine species, and biogeographic disjunctions in species assemblages. Our results indicate that temperature-mediated selection represents an early stage of marine ecological speciation in coastal regions that lack physical dispersal barriers.

Evidence of the Complexity of Gene Expression Analysis in Fish Wild Populations.

The present work examines the induction of the band 3 anion transport protein, mitogen-activated protein kinase, and lactate dehydrogenase, respectively related to osmolyte transport, cell volume regulation, and energy production in the gills of two tilapia strains exposed to either freshwater or hypersaline water. Overall, genes showed similar expression patterns between strains. However, a wild population survey across a range of natural habitats and salinities did not reveal the expected patterns. Although significant, the correlations between gene expression and salinity were slightly ambiguous and did not show any link with phenotypic differences in life history traits previously reported between the same populations. The differential expression was also not associated with the population genetic structure inferred from neutral markers. The results suggest that the differential expression observed is not the result of evolutionary forces such as genetic drift or adaptation by natural selection. Instead, it can be speculated that genes responded to various abiotic and biotic stressors, including factors intrinsic to animals. This study provides clear evidence of the complexity of gene expression analysis in wild populations and shows that more attention needs to be paid when selecting candidates as potential biomarkers for monitoring adaptive responses to a specific environmental perturbation.

The transcriptome of metamorphosing flatfish.

BACKGROUND: Flatfish metamorphosis denotes the extraordinary transformation of a symmetric pelagic larva into an asymmetric benthic juvenile. Metamorphosis in vertebrates is driven by thyroid hormones (THs), but how they orchestrate the cellular, morphological and functional modifications associated with maturation to juvenile/adult states in flatfish is an enigma. Since THs act via thyroid receptors that are ligand activated transcription factors, we hypothesized that the maturation of tissues during metamorphosis should be preceded by significant modifications in the transcriptome. Targeting the unique metamorphosis of flatfish and taking advantage of the large size of Atlantic halibut (Hippoglossus hippoglossus) larvae, we determined the molecular basis of TH action using RNA sequencing. RESULTS: De novo assembly of sequences for larval head, skin and gastrointestinal tract (GI-tract) yielded 90,676, 65,530 and 38,426 contigs, respectively. More than 57 % of the assembled sequences were successfully annotated using a multi-step Blast approach. A unique set of biological processes and candidate genes were identified specifically associated with changes in morphology and function of the head, skin and GI-tract. Transcriptome dynamics during metamorphosis were mapped with SOLiD sequencing of whole larvae and revealed greater than 8,000 differentially expressed (DE) genes significantly (p < 0.05) up- or down-regulated in comparison with the juvenile stage. Candidate transcripts quantified by SOLiD and qPCR analysis were significantly (r = 0.843; p < 0.05) correlated. The majority (98 %) of DE genes during metamorphosis were not TH-responsive. TH-responsive transcripts clustered into 6 groups based on their expression pattern during metamorphosis and the majority of the 145 DE TH-responsive genes were down-regulated. CONCLUSIONS: A transcriptome resource has been generated for metamorphosing Atlantic halibut and over 8,000 DE transcripts per stage were identified. Unique sets of biological processes and candidate genes were associated with changes in the head, skin and GI-tract during metamorphosis. A small proportion of DE transcripts were TH-responsive, suggesting that they trigger gene networks, signalling cascades and transcription factors, leading to the overt changes in tissue occurring during metamorphosis.

Diversification and coevolution of the ghrelin/growth hormone secretagogue receptor system in vertebrates.

The gut hormone ghrelin is involved in numerous metabolic functions, such as the stimulation of growth hormone secretion, gastric motility, and food intake. Ghrelin is modified by ghrelin O-acyltransferase (GOAT) or membrane-bound O-acyltransferase domain-containing 4 (MBOAT4) enabling action through the growth hormone secretagogue receptors (GHS-R). During the course of evolution, initially strong ligand/receptor specificities can be disrupted by genomic changes, potentially modifying physiological roles of the ligand/receptor system. Here, we investigated the coevolution of ghrelin, GOAT, and GHS-R in vertebrates. We combined similarity search, conserved synteny analyses, phylogenetic reconstructions, and protein structure comparisons to reconstruct the evolutionary history of the ghrelin system. Ghrelin remained a single-gene locus in all vertebrate species, and accordingly, a single GHS-R isoform was identified in all tetrapods. Similar patterns of the nonsynonymous (dN) and synonymous (dS) ratio (dN/dS) in the vertebrate lineage strongly suggest coevolution of the ghrelin and GHS-R genes, supporting specific functional interactions and common physiological pathways. The selection profiles do not allow confirmation as to whether ghrelin binds specifically to GOAT, but the ghrelin dN/dS patterns are more similar to those of GOAT compared to MBOAT1 and MBOAT2 isoforms. Four GHS-R isoforms were identified in teleost genomes. This diversification of GHS-R resulted from successive rounds of duplications, some of which remained specific to the teleost lineage. Coevolution signals are lost in teleosts, presumably due to the diversification of GHS-R but not the ghrelin gene. The identification of the GHS-R diversity in teleosts provides a molecular basis for comparative studies on ghrelin's physiological roles and regulation, while the comparative sequence and structure analyses will assist translational medicine to determine structure-function relationships of the ghrelin/GHS-R system.

Evolutionary significance and diversification of the phosphoglucose isomerase genes in vertebrates.

BACKGROUND: Phosphoglucose isomerase (PGI) genes are important multifunctional proteins whose evolution has, until now, not been well elucidated because of the limited number of completely sequenced genomes. Although the multifunctionality of this gene family has been considered as an original and innate characteristic, PGI genes may have acquired novel functions through changes in coding sequences and exon/intron structure, which are known to lead to functional divergence after gene duplication. A whole-genome comparative approach was used to estimate the rates of molecular evolution of this protein family. RESULTS: The results confirm the presence of two isoforms in teleost fishes and only one variant in all other vertebrates. Phylogenetic reconstructions grouped the PGI genes into five main groups: lungfishes/coelacanth/cartilaginous fishes, teleost fishes, amphibians, reptiles/birds and mammals, with the teleost group being subdivided into two subclades comprising PGI1 and PGI2. This PGI partitioning into groups is consistent with the synteny and molecular evolution results based on the estimation of the ratios of nonsynonymous to synonymous changes (Ka/Ks) and divergence rates between both PGI paralogs and orthologs. Teleost PGI2 shares more similarity with the variant found in all other vertebrates, suggesting that it has less evolved than PGI1 relative to the PGI of common vertebrate ancestor. CONCLUSIONS: The diversification of PGI genes into PGI1 and PGI2 is consistent with a teleost-specific duplication before the radiation of this lineage, and after its split from the other infraclasses of ray-finned fishes. The low average Ka/Ks ratios within teleost and mammalian lineages suggest that both PGI1 and PGI2 are functionally constrained by purifying selection and may, therefore, have the same functions. By contrast, the high average Ka/Ks ratios and divergence rates within reptiles and birds indicate that PGI may be involved in different functions. The synteny analyses show that the genomic region harbouring PGI genes has independently undergone genomic rearrangements in mammals versus the reptile/bird lineage in particular, which may have contributed to the actual functional diversification of this gene family.

European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation.

The European sea bass (Dicentrarchus labrax) is a temperate-zone euryhaline teleost of prime importance for aquaculture and fisheries. This species is subdivided into two naturally hybridizing lineages, one inhabiting the north-eastern Atlantic Ocean and the other the Mediterranean and Black seas. Here we provide a high-quality chromosome-scale assembly of its genome that shows a high degree of synteny with the more highly derived teleosts. We find expansions of gene families specifically associated with ion and water regulation, highlighting adaptation to variation in salinity. We further generate a genome-wide variation map through RAD-sequencing of Atlantic and Mediterranean populations. We show that variation in local recombination rates strongly influences the genomic landscape of diversity within and differentiation between lineages. Comparing predictions of alternative demographic models to the joint allele-frequency spectrum indicates that genomic islands of differentiation between sea bass lineages were generated by varying rates of introgression across the genome following a period of geographical isolation.

The impact of environmental degradation on reproduction of the black-chinned tilapia Sarotherodon melanotheron from various coastal marine, estuarine and freshwater habitats.

There is growing evidence that climate change has greatly altered environmental conditions in many aquatic ecosystems over the last decades, leading to changes in fish distribution and life history traits. Recent works conducted in Senegalese and Gambian coastal marine, estuarine and freshwater ecosystems have shown important changes in the intensity, frequency and breeding timing of wild population of Sarotherodon melanotheron in response to changes in salinity regimes. In addition to salinity, this study investigates the potential influences of other environmental factors that have received less attention on the reproduction of S. melanotheron from three different aquatic ecosystems in Senegal. The results demonstrate that day-length and temperature affect sexual maturity in both males and females from Guiers Lake and Hann Bay, but no such effects were found in upstream of the Saloum Estuary, where the spawning activity seems to be under the synergetic control of rainfall and salinity that apparently predominate over all the other ambient factors. This study demonstrated for the first time that aside from photoperiod and temperature, rainfall also influences fish reproduction in Hann Bay probably through its effects on water quality. Furthermore, our results confirm previous findings that changes in salinity regimes resulting from seasonal variations in precipitation and evaporation are the ultimate causes of temporal changes in the reproductive activity of the species in this estuary.

Genomic characterization of the European sea bass Dicentrarchus labrax reveals the presence of a novel uncoupling protein (UCP) gene family member in the teleost fish lineage.

BACKGROUND: Uncoupling proteins (UCP) are evolutionary conserved mitochondrial carriers that control energy metabolism and therefore play important roles in several physiological processes such as thermogenesis, regulation of reactive oxygen species (ROS), growth control, lipid metabolism and regulation of insulin secretion. Despite their importance in various physiological processes, their molecular function remains controversial. The evolution and phylogenetic distribution may assist to identify their general biological function and structure-function relationships. The exact number of uncoupling protein genes in the fish genome and their evolution is unresolved. RESULTS: Here we report the first characterisation of UCP gene family members in sea bass, Dicentrarchus labrax, and then retrace the evolution of the protein family in vertebrates. Four UCP genes that are shared by five other fish species were identified in sea bass genome. Phylogenetic reconstitution among vertebrate species and synteny analysis revealed that UCP1, UCP2 and UCP3 evolved from duplication events that occurred in the common ancestor of vertebrates, whereas the novel fourth UCP originated specifically in the teleost lineage. Functional divergence analysis among teleost species revealed specific amino acid positions that have been subjected to altered functional constraints after duplications. CONCLUSIONS: This work provides the first unambiguous evidence for the presence of a fourth UCP gene in teleost fish genome and brings new insights into the evolutionary history of the gene family. Our results suggest functional divergence among paralogues which might result from long-term and differential selective pressures, and therefore, provide the indication that UCP genes may have diverse physiological functions in teleost fishes. Further experimental analysis of the critical amino acids identified here may provide valuable information on the physiological functions of UCP genes.

Comparative analysis of reproductive traits in black-chinned tilapia females from various coastal marine, estuarine and freshwater ecosystems.

The black-chinned tilapia Sarotherodon melanotheron is a marine teleost characterised by an extreme euryhalinity. However, beyond a certain threshold at very high salinity, the species exhibits impaired growth and precocious reproduction. In this study, the relationships between reproductive parameters, environmental salinity and condition factor were investigated in wild populations of this species that were sampled in two consecutive years (2003 and 2004) from three locations in Senegal with different salinities: Guiers lake (freshwater, 0 psu), Hann bay (seawater, 37 psu) and Saloum estuary (hypersaline water, 66-127 psu). The highest absolute fecundity and spawning weight were recorded in seawater by comparison to either freshwater or hypersaline water whereas the poorest condition factors were observed in the most saline sampling site. These results reflect higher resource allocation to the reproduction due to the lowest costs of adaptation to salinity in seawater (the natural environment of this species) rather than differences in food resources at sites and/or efficiency at foraging and prey availability. Fecundities, oocyte size as well as spawning weight were consistent from year to year. However, the relative fecundity in the Saloum estuary varied significantly between the dry and rainy raisons with higher values in the wet season, which seems to reflect seasonal variations in environmental salinity. Such a reproductive tactic of producing large amounts of eggs in the rainy season when the salinity in the estuary was lower, would give the fry a better chance at survival and therefore assures a high larval recruitment. An inverse correlation was found between relative fecundity and oocyte size at the two extreme salinity locations, indicating that S. melanotheron has different reproductive strategies in these ecosystems. The adaptive significance of these two reproductive modes is discussed in regard to the heavy osmotic constraint imposed by extreme salinities and high inter-specific competition.

Directed sequencing and annotation of three Dicentrarchus labrax L. chromosomes by applying Sanger- and pyrosequencing technologies on pooled DNA of comparatively mapped BAC clones.

Dicentrarchus labrax is one of the major marine aquaculture species in the European Union. In this study, we have developed a directed-sequencing strategy to sequence three sea bass chromosomes and compared results with other teleosts. Three BAC DNA pools were created from sea bass BAC clones that mapped to stickleback chromosomes/groups V, XVII and XXI. The pools were sequenced to 17-39x coverage by pyrosequencing. Data assembly was supported by Sanger reads and mate pair data and resulted in superscaffolds of 13.2 Mb, 17.5 Mb and 13.7 Mb respectively. Annotation features of the superscaffolds include 1477 genes. We analyzed size change of exon, intron and intergenic sequence between teleost species and deduced a simple model for the evolution of genome composition in teleost lineage. Combination of second generation sequencing technologies, Sanger sequencing and genome partitioning strategies allows "high-quality draft assemblies" of chromosome-sized superscaffolds, which are crucial for the prediction and annotation of complete genes.

Comparative analysis of intronless genes in teleost fish genomes: insights into their evolution and molecular function.

This study assessed the relationship between the occurrence and function of intronless or single exon genes (SEG) in the genome of five teleost species and their phylogenetic distance. The results revealed that Takifugu rubripes, Tetraodon nigroviridis, Oryzias latipes, Gasterosteus aculeatus and Danio rerio genomes are respectively comprised of 2.83%, 3.42%, 4.49%, 4.35% and 4.02% SEGs. These SEGs encode for a variety of family proteins including claudins, olfactory receptors and histones that are essential for various biological functions. Subsequently, we predicted and annotated SEGs in three European sea bass, Dicentrarchus labrax chromosomes that we have sequenced, and compared results with those of stickleback (G. aculeatus) homologous chromosomes. While the annotation features of three D. labrax chromosomes revealed 78 (5.30%) intronless genes, comparisons with G. aculeatus showed that SEG composition and their order varied significantly among corresponding chromosomes, even for those with nearly complete synteny. More than half of SEGs identified in most of the species have at least one ortholog multiple exon gene in the same genome, which provides insight to their possible origin by retrotransposition. In spite of the fact that they belong to the same lineage, the fraction of predicted SEGs varied significantly between the genomes analyzed, and only a low fraction of proteins (4.1%) is conserved between all five species. Furthermore, the inter-specific distribution of SEGs as well as the functional categories shared by species did not reflect their phylogenetic relationships. These results indicate that new SEGs are continuously and independently generated after species divergence over evolutionary time as evidenced by the phylogenetic results of single exon claudins genes. Although the origin of SEGs cannot be inferred directly from the phylogeny, our results provide strong support for the idea that retrotransposition followed by tandem duplications is the most probable event that can explain the expansion of SEGs in eukaryotic organisms.

A Comparative BAC map for the gilthead sea bream (Sparus aurata L.).

This study presents the first comparative BAC map of the gilthead sea bream (Sparus aurata), a highly valuated marine aquaculture fish species in the Mediterranean. High-throughput end sequencing of a BAC library yielded 92,468 reads (60.6 Mbp). Comparative mapping was achieved by anchoring BAC end sequences to the three-spined stickleback (Gasterosteus aculeatus) genome. BACs that were consistently ordered along the stickleback chromosomes accounted for 14,265 clones. A fraction of 5,249 BACs constituted a minimal tiling path that covers 73.5% of the stickleback chromosomes and 70.2% of the genes that have been annotated. The N50 size of 1,485 "BACtigs" consisting of redundant BACs is 337,253 bp. The largest BACtig covers 2.15 Mbp in the stickleback genome. According to the insert size distribution of mapped BACs the sea bream genome is 1.71-fold larger than the stickleback genome. These results represent a valuable tool to researchers in the field and may support future projects to elucidate the whole sea bream genome.

Analysis of single nucleotide polymorphisms in three chromosomes of European sea bass Dicentrarchus labrax.

Single nucleotide polymorphisms (SNPs) are believed to contain relevant information and have been therefore extensively used as genetic markers in population and conservation genetics, and molecular ecology studies. This study reports on the identification of potential SNPs in a diploid European sea bass Dicentrarchus labrax genome by using reference sequences from three assembled chromosomes and mapping all WGS datasets onto them (3× Sanger, 3× 454 and 20× SOLEXA). A total of 20,779 SNPs were identified over the 1469 gene loci and intergenic space analysed. Within chromosomes the occurrence of SNPs was the lowest in exons and higher in introns and intergenic regions, which may be explained by the fact, that coding regions are under strong selective pressure to maintain their biological function. The ratio of nonsynonymous to synonymous mutations was smaller than one for all the chromosomes, suggesting that most of deleterious nonsynonymous mutations were eliminated by negative selection. SNPs were not uniformly distributed over the chromosomes. Two chromosomes exhibited large regions with extremely low SNP density, which might represent homozygous regions in the diploid genome. The results of this study show how SNP detection can take profit from sequencing a single diploid individual, but also uncover the limits of such an approach. SNPs that have been identified will support marker development for genetic linkage mapping, population genetics and aquaculture related questions in general.

Differential expression of the heat shock protein Hsp70 in natural populations of the tilapia, Sarotherodon melanotheron, acclimatised to a range of environmental salinities.

BACKGROUND: The relationship between environmental variation and induction of heat shock proteins (Hsps) has been much documented under experimental conditions. However, very little is known about such induction in natural populations acclimatised to prevailing environmental conditions. Furthermore, while induction of stress proteins has been well documented in response to environmental contaminants and thermal stressors, little is known about whether factors, such as extreme salinity, are also potential inductors. The black-chinned tilapia Sarotherodon melanotheron is unusual for its ability to colonise estuarine environments in West Africa that are characterised by extremely high salinities. The relationships between mRNA levels of the 70 kDa heat shock protein (Hsp70) and Na+, K+-ATPase1alpha (Naka) in the gills, environmental salinity, and a life-history trait (condition factor) were investigated in wild populations of this species sampled from three locations in the Saloum estuary, at salinities ranging from 40 to 100 psu. RESULTS: The highest Hsp70 and Naka mRNA levels, and the poorest condition factors were recorded in the most saline sampling site (100 psu). The Hsp70 and Naka mRNA were correlated amongst themselves and showed a direct positive correlation with environmental salinity, and a negative correlation with fish condition factor. Thus, the Hsp70 is constitutively overexpressed by S. melanotheron acclimatised to extreme hypersalinity. CONCLUSIONS: These results indicate that, although S. melanotheron can colonise extremely saline environments, the overexpression of Hsp70 combined with the higher Naka mRNA expression reveals that this represents a chronic stress. The induction of Hsp70 was, therefore, a biomarker of chronic hyper-osmotic stress which presumably can be linked to the impaired growth performance and precocious reproduction that have been demonstrated in the populations at the extremely saline sites.

Transcriptional responses of the black-chinned tilapia Sarotherodon melanotheron to salinity extremes.

Sarotherodon melanotheron is one of the most euryhaline teleosts able to withstand variations in environmental salinity ranging from freshwater (FW) to 130‰ hyper-saline waters (HSW). Although significant progress has been made in exploring the cellular and molecular changes that accompany salinity adaptation in teleosts, little is known about the effects of long-term acclimation to HSW. We sought to identify in this tilapia species the genes whose transcription is induced by long-term acclimation either to HSW or FW. Two subtractive cDNA libraries were made from gills of fish acclimated for 45 days to either condition, with 320 partial cDNA sequences encoding proteins potentially involved in the response to the two salinity extremes. The ESTs comparisons with genomic databases allowed putative functions to be attributed to 197 of these genes. The suppression subtractive hybridisation (SSH) results were validated by Real-time PCR for 13 candidate genes having presumably a role in osmoregulation, supplemented by Na(+), K(+)-ATPase α-subunit and carbonic anhydrase, two genes known to be implicated in this function. In fish acclimated to both salinity extremes, the functional category of cellular process was the predominant one, which may indicate high cellular turnover rates in FW and HSW-adapted fish. The acclimation to FW and HSW also appeared to trigger the expression of genes involved in transport activity, biological regulation and metabolic processes, at a higher level in fish acclimated to HSW, suggesting higher metabolic activity in this situation. These results are a first step towards the identification of key molecular processes involved in the fish acclimation to extreme salinities.